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bJ!) Abstract 

Callan, Giddings, Harvey and Stroininger (CGHS) previously in- 

^ troduced a two-dimensional semiclassical model of gravity coupled to 

\^ a dilaton and to matter fields. Their model yields a system of field 

zj^ equations which may describe the formation of a black hole in gravita- 

*y-v tional collapse as well as its subsequent evaporation. Here we present 

1^ an approximate analytical solution to the semiclassical CGHS field 

^^ equations. This solution is constructed using the recently-introduced 

^^ formalism of flux-conserving hyperbolic systems. We also explore the 

. . asymptotic behavior at the horizon of the evaporating black hole. 

> 

H 1 Introduction 

The semiclassical theory of gravity treats spacetime geometry at the classi- 
cal level but allows quantum treatment of the various fields which reside in 
spacetime. This theory asserts that a quantum field living on a black-hole 
(BH) background will usually be endowed with non-trivial fluxes of energy- 
momentum. These fluxes, represented by the renormalized stress-Energy ten- 
sor Taf^, originate from the field's quantum fluctuations, and typically they 
do not vanish even in the (incoming) vacuum state. The Hawking radiation 



[1], and the consequent black-hole (BH) evaporation, are perhaps the most 
dramatic manifestations of these quantum fluxes. 

In the framework of semiclassical gravity the spacetime reacts to the quan- 
tum fluxes via the Einstein equations, which now receive the extra quantum 
contribution T^^ at their right-hand side. The mutual interaction between ge- 
ometry and quantum fields thus takes its usual General- Relativistic schematic 
form: The renormalized stress-Energy tensor Tq/j (say in the vacuum state) 
is dictated by the background geometry, and the latter is affected by Tq,^ 
through the Einstein equations. In principle, this evolution scheme allows 
systematic investigation of the spacetime of an evaporating BH. 

It turns out, however, that the calculation of Tq^(x) for a prescribed back- 
ground metric Qafsix) is an extremely hard task in four dimensions. Never- 
theless, in two-dimensional (2D) gravity the situation is remarkably simpler. 
There are two energy-momentum conservation equations, and the trace T^ 
is also known (the "trace- anomaly"). These three pieces of information are 
just enough for determining the three unknown components of T^/J- It is 
therefore possible to implement the evolution scheme outlined above in 2D- 
gravity, and to formulate a closed system of field equations which describe 
the combined evolution of both spacetime and quantum fields. 

Callan, Giddings, Harvey and Strominger (CGHS) [2J introduced a for- 
malism of 2D gravity in which the metric is coupled to a dilaton field (j) 
and to a large number A^ of identical massless scalar fields. They added to 
the classical action an effective term oc A^ which gives rise to the semiclassi- 
cal trace-anomaly contribution, and thereby automatically incorporates the 
renormalized stress-Energy tensor T^^ into spacetime dynamics. The evo- 
lution of spacetime and fields is then described by a closed system of field 
equations. They considered the scenario in which a BH forms in the grav- 
itational collapse of a thin massive shell, and then evaporates by emitting 
Hawking radiation. Their main goal was to reveal the end-state of the evap- 
oration process, in order to address the information puzzle. 

The general analytical solution to the CGHS field equations is not known. 
Nevertheless, these equations can be explored analytically [2] as well as nu- 
merically j3l El E]. The global structure which emerges from these studies 
is depicted in Fig. 1: The shell collapse leads to the formation of a BH. A 
spacelike singularity forms inside the BH, at a certain critical value of the 
dilaton [3]. (The BH interior may be identified as the set of all events from 
which all future-directed causal curves hit the singularity.) The BH interior 



and exterior are separated by an outgoing null ray, which serves as an event 
horizon. Also an apparent horizon forms along a timelike line outside the BH 
(it is characterized by a local minimum oi R = e"^*^ along outgoing null rays). 
The event and apparent horizons are denoted by "EH" and "AH" in Fig. 1. 
Both horizons steadily "shrink" it time (namely R decreases), exhibiting the 
BH evaporation process. At a certain point the apparent horizon intersects 
with the spacelike singularity (and with the event horizon). This intersection 
event (denoted "P" in Fig. 1) appears to be a naked singularity, visible to 
far asymptotic observers. It may be regarded as the "end of evaporation" 
point. 

The spacelike singularity which develops inside a CGHS evaporating BH 
was first noticed by Russo, Susskind, and Thorlacius [3j. Its local structure 
was recently studied in some detail [7], by employing the homogeneous ap- 
proximation. This singularity marks the boundary of predictability of the 
semiclassical CGHS formalism. About two years ago Ashtekar, Taveras and 
Varadarajan [S] proposed a quantized version of the CGHS model, in which 
the dilaton and metric are elevated to quantum operators. This quantum 
formulation of the problem seems to resolve the semiclassical singularity in- 
side the BH [8j, and thereby to shed new light on the information puzzle. In 
a very recent paper [9] we applied a more simplified, "Minisuper space-like", 
quantum approach to the spacetime singularity inside a CGHS evaporating 
BH. We obtained a (bounce-like) extension of semiclassical spacetime beyond 
the singularity. 

The goal of this paper is to construct an approximate analytical solution 
to the semiclassical CGHS field equations, which will satisfactorily describe 
the BH formation and evaporation. This approximate solution applies as long 
as the BH is macroscopic, and as long as the spacetime region in consideration 
is "macroscopic" too — namely, R is sufficiently large (i.e. not too close to 
the singularity). In such macroscopic regions, the semiclassical effects are 
weak in a local sense. Our solution is "first-order accurate" in this local 
sense; Namely, the local errors in the second-order field-equation operators 
(applied to the approximate solution) are quadratic in the magnitude of local 
semiclassical effects (this is further discussed in Sec. ^. 

The method we use here for constructing our approximate solution is 
based on the formalism oi flux- conserving systems jlO]. This formalism deals 
with a special class of semi-linear second-order hyperbolic systems in two di- 
mensions. This class was first introduced in Ref. [TT], and was subsequently 
explored more systematically in Ref. jl2] (see also [13]). In particular, a 



flux- conserving hyperbolic system admits a rich family of single-flux solu- 
tions, to which we refer as Vaidya-like solutions (see Sec. [S]), and for which 
the field equations reduce to an ordinary differential equation (ODE). Very 
recently we demonstrated [lOj that after transforming to new field variables, 
the semiclassical CGHS equations take a form which is approximately flux- 
conserving. Here we take advantage of this property and use the formalism 
of flux-conserving systems to construct the approximate analytic solution to 
the CGHS field equations. 

Our approximate solution is presented in Sec. |6| It involves a single 
function denoted by H, which is defined through a certain ODE. In order 
to make practical use of this approximate solution, one must have at his 
disposal the solution for this ODE. In Appendix[X]we provide an approximate 
analytical solution to this ODE. 

In Sec. |2] we briefly summarize the CGHS model, expressing the semi- 
classical field equations in convenient variables R, S (already used in Refs. 
[TTllTtlTU]). In particular we discuss the classical model of a collapsing shell — 
which in fact provides the initial conditions to the semiclassical problem of 
BH evaporation. In Sec. [slwe introduce the new field variables {W, Z), which 
turn the field equations into a more standard form (with no first-order deriva- 
tives). Then we observe the approximate flux-conserving system obtained in 
these new variables at the large- i? limit. In Sec. |4]we construct the ingoing 
Vaidya-like solution (a "single-flux" ingoing solution) of the flux- conserving 
system, which constitutes the core to our approximate solution. This core 
solution is complemented in Sec. |5] by adding to it a weak outgoing com- 
ponent. The extra outgoing component must be added in order to correctly 
satisfy the initial conditions at the collapsing shell (and it is this compo- 
nent which eventually gives rise to the outgoing Hawking radiation). In Sec. 
|6] we summarize our approximate solution, and also present an alternative 
approximate expression for S. In Sec. [7] we introduce two useful gauges: 
the " shifted-Kruskal " gauge, and the semiclassical Eddington gauge. The 
validity of the constructed solution is verified in Sec. [8j We first describe 
the magnitude of the local error in the field equations, which turns out to 
be quadratic in the local magnitude of semiclassical effects (assumed to be a 
small quantity). Subsequently we verify that the initial data at the collapsing 
shell are precisely matched by our approximate solution, and the same for 
the initial data at past-null-infinity (PNI). 

In Sec. [9] we discuss the horizon of the semiclassical BH and a few related 
issues. The horizon (or "event-horizon") is the outgoing null line which 



separates the BH interior and exterior. We present the behavior of R and 
S along the event horizon and also compute the influx T„j, into the BH. We 
also discuss the location of the apparent horizon (deflned by a mininium of R 



along an outgoing null ray as mentioned above). Then in Sec. 10 we analyze 
the asymptotic behavior at PNI, and also remark briefly on the asymptotic 
behavior at future null inflnity (FNI). Finally, in Sec. 11 we summarize our 
main results. 

2 Background: The CGHS model 

2.1 Action and field equations 

The CGHS model P] consists of 2D gravity coupled to a dilaton and to a 
large number A^ ^ 1 of identical (free, minimally-coupled, massless) scalar 
flelds fi. Throughout this paper we express the metric in double- null coor- 
dinates u, V for convenience, namely 

ds^ = -e^Pdudv. (1) 

The action then takes the form 



1 /• r 1 ^ 



(2) 



where K = N/12. The term A^ denotes a cosmological constant. We shall set 
A = 1 throughout. This is achieved by a change of variable p ^ p' = p+ln(A), 
which annihilates A but does not affect the fleld equations otherwise. This 
setting actually amounts to a choice of basic length unit. 
The scalar flelds all satisfy the trivial fleld equation 

The CGHS scenario consists of an imploding thin massive shell of mass Mq 
which moves along an ingoing null line v = const = vq, forming a black 
hole. This shell is composed of (f-derivatives of) the scalar flelds f], which 
are concentrated as a (5-f unction-like distribution at v = Wq. The solution at 
f < f is the trivial, flat, vacuum solution (see below). The main objective 
of this paper is the semiclassical solution which takes place at f > wq. By 



assumption, in this region too no incoming scalar waves are present. There- 
fore the solution of Eq. ([S]) throughout the relevant domain v > Vq (as well 

as at t> < vq) is 

Mu,v) = 0. (4) 

The remaining field equations consist of two evolution equations for the 
fields 0, p, as well as two constraint equations (see |2]). To bring these equa- 
tions to a simpler form we define new field variables 

R^e''^, S^2{p -<!>), (5) 

following Ref. [TT]. The field equations for R and S then take the form 

Ruv = -e^ - Kp^uv, (6) 

S^uv = Kp^uv/R, (7) 

where 

p={S-\nR)/2 (8) 

is to be substituted. The constraint equations become (after substituting 
ft = for the scalar-fields) 

R,VV R,V^,V 'T -LyV ^ ^t {^^) 

where Tuu, T^y are the semiclassical fluxes in the two null directions, given by 

Tuu = K[p^uu- P^^u + ^uiu)], (11) 

f;, = K [p,„, - pI + z.iv)] . (12) 

The functions z„(m), z^{v) are initial functions which encode the information 
about the system's quantum state. Following CGHS we consider here an in- 
coming vacuum state at PNI (apart from the imploding massive shell, which 
is presumably encoded in the fields fi already at the classical level). Corre- 
spondingly the functions z^iu), z^{v) are determined by the requirement that 
Tuu and T^u vanish at PNI. 



^In a 2D spacetime there are two sectors of PNI, a "right PNI" and a "left PNI". z„(m) 
is defined by the demand T^u = at left PNI, and z^{v) by demanding Ty^ = at right 
PNI. Similarly there are two sectors of FNI, a right one and a left one. Throughout this 
paper, by "PNI" and "FNI" we shall always refer to the right sectors of PNI and FNI 
(unless stated otherwise), as is also illustrated in Fig. 1. 
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It will sometimes be useful to re-express the system of evolution equations 
(GIT) in its standard form, in which R^uv and S^uv are explicitly given in terms 



of lower-order derivatives: 

■"'' ~ '"" 2{R-K) ~ '"" '''2R{R-Ky ^'^^> 

This form makes it obvious that the evolution equations become singular 
when R = K, and also at i? = 0. This singularity was studied in some detail 
inRef. [7]. 

Gauge freedom 

In a coordinate transformation u — )■ u'{u),v — )■ v'{v) the dilaton scalar field 
(p is unchanged, but p changes as 

, 1 / du' , dv'\ .^^, 

(as may be deduces from the coordinate transformation of the metric com- 
ponent guv = ~(l/2)e^''). From the definition (pi) of i? and S it is obvious 
that i? is a scalar, and S changes in a coordinate transformation like 2p: 

R' = R, 5' = 5-ln^-ln^. (16) 

du dv 

Below we shall often provide expressions for 5* in certain specific gauges. 
In such cases we shall use the notation S'[...], with the specific u, v coordinates 
specified in the squared brackets. The same notation will apply to p (and 
also to the gauge-dependent quantity Z introduced in the next section). 

2.2 Classical solutions 

The classical solutions are obtained by setting K = 0, leading to T^u = T^^ = 
0. The vacuum field equations then reduce to the evolution equations 

R,uv ^ 6 ; ^,uv ^ U (-'-' j 



and the constraint equations 

The general solution of these equations may be easily constructed. It takes 
the form 

Riu, v) = M- Ru{u)R,{v) , S{u, v) = \n{Ru,uRv,v), (19) 

where M is an arbitrary constant, and Ru{u) and Rv{v) are any monotonically- 
increasing functions of their arguments (with non- vanishing derivatives). The 
classical solution may thus look at first glance as a rich class depending on 
two arbitrary functions. However, these two functions merely reflect the 
gauge freedom. To fix this freedom we may use the Kruskal-like coordinates 
U = Ru{u),V = Rv{v), after which the general solution takes the simple 
explicit form 

R = M-UV , Spy] = 0. (20) 

The sub-index " [U, V] " , recall, indicates that this expression for 5* only ap- 
plies in a specific gauge, the one associated with the Kruskal U, V coordinates. 



The representation (20) makes it obvious that the classical vacuum solu- 
tion is a one-parameter family, parametrized by the mass M. We shall refer 
to it as the Schwarzschild-like solution. 

For M > the spacetime contains a BH, whose causal structure resembles 
that of the four- dimensional Schwarzschild spacetime (see [2]). The event 
horizon is located at f/ = 0, and the past (or "white-hole") horizon at V^ = 0. 
Inside the BH {U, V" > 0) there is a spacelike R = singularity (0, p diverge) 
at UV = M. For negative M there is a naked, timelike, R = singularity 
instead of a BH. (The M = case is considered below.) 

Another useful gauge is the Eddington-like gauge, obtained by the trans- 
formation G 

Ue = -\n{-U), Ve = \TiV (f/<0,1^>0). (21) 

The solution then takes the form 

i? = M + e^^-"% 5[„,,.,] = t;e - «e. (22) 



^The coordinate u^, introduced here is the classical outgoing Eddington coordinate. It 
should not be confused with the seniiclassical outgoing Eddington coordinate u introduced 
in Sec. It] On the other hand the ingoing Eddington coordinate Ve (denoted later by v) is 
common to the classical and semiclassical solutions. 



Note that these coordinates only cover the BH exterior. Analogous Eddington- 
like coordinates (with appropriate sign changes) may also be constructed for 
the BH interior, but we shall not need these internal coordinates here. 

The special case M = yields the Minkowski- like solution. In Kruskal 
coordinates it takes the form 

R=-UV, Siu,v]=0. (23) 

In Eddington coordinates one obtains R = e^'=~"'=, S[ue,ve] = Ve — u^ and hence 
P[ue,v^] = 0. This demonstrates that the M = solution is flat, and also 
indicates that the Eddington-like coordinates {ue,Ve) correspond in this case 
to the standard, flat, null coordinates of 2D Minkowski spacetime. 



In the general case M 7^ one finds from Eqs. (22) and l\8 



and the associated metric ds'^ = —e'^^dudv yields non-vanishing curvature. 
However, at large Ve — Ue this expression reduces to 

P[ueM = - 2^ve-u, ~ ~2R' ^ > 

which vanishes as fg — Ug — ?■ 00. That is, in Eddington coordinates p vanishes 
at (right) PNI {ue -^ —00), FNI {ve — ?■ 00), and spacelike infinity {ue — >■ 
— oo,fe — i- 00). The Schwarzschild-like solution is thus asymptotically-flat, 
and {ueiVe) serve as asymptotically-flat coordinates in this spacetime. In 
particular, v^ and Uf. respectively coincide with the affine parameter along 
PNI and FNI. 

Classical collapsing shell 

The physical scenario which concerns us here is the formation of a BH by 
the collapse of a thin shell, carrying a mass Mq > 0, which propagates along 
an ingoing null line v = Vq. The classical solution describing this scenario 
is obtained by continuously matching a Minkowski-like solution sX v < Vq 
with a Schwarzschild-like solution at f > i;o- To describe the matching we 



start with the Kruskal-gauge solution (20 ), along with its M = analog (23 ) 



Obviously the i?-functions in these two expressions do not properly match 
(at any V). To allow for continuous matching, we shift the U coordinate in 
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the Minkowski-like domain f < f o , defining U = U + AU, where AU is a 
constant. Note that S is unchanged in this coordinate transformation. We 
obtain the Minkowski-like solution in its new form: 

R={AU-U)V, %y]=0 {v<vo). (26) 



This solution can now be matched to Eq. (20), by equating VAU to Mq at 
the null shell. Omitting the hat from U and re-naming this new Minkowski- 
like Kruskal coordinate by U, we obtain the matched solution 

S„„-, = 0, mf) = ('^^ -''''' ^^^^°' (27) 

^ ' ' ^ ^ [Mo -UV V>Vo, ^ 

where Vq is the shell's \^- value (we assume Vq > 0), and continuity implies 
AU = Mq/Vq. 

To formulate the semiclassical initial data (next subsection) we find it 
most convenient to work with Eddington v and Kruskal U . The collapsing- 
shell solution then takes the form 






where f^ = In Vq. 

Fixing the t>-gauge 

Throughout the semiclassical analysis below, we find the Eddington Ve to be 
the most convenient w-coordinate, and no other u-coordinate will be used. 
To simplify the notation, in the rest of this paper we shall simply denote Vg 
by f. 

For the w-gauge we shall occasionally use several ^-coordinates at various 
stages of the construction. The value of S (like p and Z) in specific w-gauges 
will be denoted by S*!..,] — with the specific u-coordinate only specified in the 
squared brackets — and it will always refer to the Eddington v coordinate. No 
confusion should arise, because no f-coordinate other than Eddington will 
be used throughout the rest of the paper. 

To practice the new notational rules we re- write the above classical collapsing 



shell solution (28) in the new notation: 



V < Vq, 



S,u,=v, R{U,v) = \[^^ "/l'^ —- (29) 

I Mo — Ue V > Vq, 
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where f q = In Vq denotes the shell's location in Eddington-f , and 



The collapsing-shell spacetime (29) is asymptotically-flat, just like the 



pure Schwarzschild-like and Minkowski-like solutions. Also, the Eddington 
coordinate v serves as an afline parameter at PNI. [This may be verifled by 
switching to the Eddington coordinate u^ = — ln(— f/) and noting that p[u^] 
vanishes at the PNI limit Ue — )■ — oo.] 

2.3 Initial data for the semiclassical solution 

Initial conditions at the Shell 

We return now to the semiclassical problem of collapsing shell. The semi- 
classical effects vanish at the portion i; < f o of the collapsing-shell spacetime, 
owing to its flatness, [j Therefore, the portion f < fo is correctly described 



by Eq. (29) even in the semiclassical problem. 

The portion v > Vq oi the collapsing-shell spacetime will be profoundly 
modifled by semiclassical effect (as expressed for example by the BH evapo- 
ration). However, by continuity, the initial conditions for R and S at v = vq 
will still be determined by matching to the same, unmodified, Minkowski-like 
solution at f < t;o, and will therefore be exactly the same as in the classi- 



cal solution (29). Using again the Kruskal-[/ coordinate (and, recall, the 



Eddington-t; coordinate as usual), these shell initial conditions take the form 

RoiU) = R{U, vo) = Mo - Ue""" (30) 

and 

Sl^\u)^S[u]iU,vo) = vo. (31) 

Initial conditions at Past null infinity 



The semiclassical evolution equations (617) also require initial conditions at 



an outgoing null ray, which is most conveniently taken to be the PNI bound- 
ary. The situation here is somewhat similar to that of the initial data at the 
shell: Owing to asymptotic flatness of the collapsing-shell spacetime [and 



•^Flatness means that the term p^y in the semiclassical evolution equations (6l7l vanish, 
so these equations reduce to the classical ones. 
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to the lack of any strong-field region at the causal past of PNI (unlike the 
situation at FNI)], the asymptotic behavior at PNI is just the classical one. 
Specifically one finds the initial conditions 

R^Mo-Ue\ S[u]=v (PNI). (32) 

in the domain v > vq. In particular, the influx T^^ should vanish at PNI 
(corresponding to an incoming vacuum state). 

2.4 The role of K in semiclassical dynamics 

The semiclassical effects originate from the last term Kp^p^v in fhe CGHS 
action ([2]). The parameter K = N/12 thus determines the overall magnitude 
of semiclassical effects. It is important to note, however, that K may be 
factored out from CGHS dynamics by a simple shift/rescaling of variables. 
In the transformation 

N^cN, K^cK, 0->0-^lnc (33) 

[with A, p and /, unchanged (and with all scalar fields fi identical)], the 
action is merely multiplied by c, hence the field equations are unaffected. 
One can therefore factor out K this way by choosing c = 1/K. 

In the relevant domain v > Vq the scalar fields fi vanish anyway, so the 

N n 

term — | XI fi,ufi,v is absent from the action (2). In terms of the new variables 

i=l ' ' |_ 

R, S, the rescaling law takes the form 

K^cK, R^cR, S^S + \nc. (34) 

This scaling makes it obvious that the relative magnitude of semiclassical 
effects in various regions of spacetime will not be determined by iT or i? 
separately, but only through (dimensionless r« combinations like K/R or 
K/e^, which are invariant to the rescaling. n 

One can easily verify that in the above scaling transformation the BH 
original mass Mq is multiplied by c. It is a common wisdom that when a BH 



''Note that after A is set to 1 [see discussion following Eq. ([2|] all model's quantities 
become dimensionless. This includes R, and also the BH mass Mq. (The latter is naturally 
linked to R — for example, through the horizon's R value.) 

'^It is important not to confuse here between two different issues related to the magni- 
tude of K: First, whether the semiclassical treatment is applicable or not. Second, within 
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is "macroscopic", semiclassical effects will be locally weak r] (except in the 
neighborhood of the singularity). The above scaling law makes it obvious 
that whether a BH may be regarded as "macroscopic" or not, would only 
depend on the ratio of Mq and K: A CGHS BH should be regarded as 
macroscopic if Mq ^ K — which we indeed assume throughout this paper. 

In principle, the above scaling allows us to set K = 1 in the field equations. 
We find it more convenient, however, to leave K in the equations untouched. 
The parameter K serves as a "flag" marking the terms of semiclassical origin 
in the various equations. Also, in the equations below K always appears 
through combinations like K/R, K/Rq, K/Mq (or sometimes with R replaced 
by the variable W ^ R introduced below). We shall assume throughout this 
paper that the BH is macroscopic (Mq ^ K), and deal with spacetime regions 
satisfying R ^ K. This will allow us to expand various expressions to first 
order in small quantities such as K/R or alike. All these expansions may 
conveniently be handled formally as expansions in K (or in the parameter 
q = K/4 introduced below) — though one should bear in mind that the small 
parameter in the expansion is not K itself, but the combinations K/R, K/Mo, 
etc. 

3 Field redefinition and flux-conserving formu- 
lation 

3.1 Field redefinition 



The evolution equations (617) may look quite simple at first glance. However 



to close the system one must substitute Eq. (^ for p, which makes the 
equations rather messy. Bringing these equations to their standard form, 
one ends up with Eqs. (13p4). In addition to their rather complicated form. 



the semiclassical formulation, K determines the magnitude of semiclassical effects. CGHS 
pointed out [2 that the semiclassical treatment is only valid if A' ^ 1. Here we address 
the second issue. Namely, we assume that the condition iiT ^ 1 is satisfied, and explore 
the scaling law which characterizes the magnitude of semiclassical effects (and its relation 
to the scale of other variables like R) . 

^Here, again, it is important to distinguish between two different aspects of "macro- 
scopicality": (i) Whether the semiclassical theory is applicable to the BH evaporation; 
and (ii) whether, within the semiclassical formulation, the semiclassical effects are locally 
weak. The discussion here pertains to the second aspect. To satisfactorily deal with (i) 
one has to further assume K ':^ \ (see also previous footnote) . 
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these second-order equations also have the inconvenient property of being 
exphcitly dependent on first-order derivatives R^u, R,v 

To get rid of this undesired dependence upon R^u and R^, we transform 
from R and S to new variables W, Z defined by |J^ 



W{R) = ^R{R - K) -K\n{y/R+VR -K)+K (^ + ln2\ 



and 
where 



Z = S + AZ{R), 



AZ{R) = — (R-W)- In R. 
K 



(35) 



(36) 



(37) 



With these new variables the evolution equations take the schematically- 
simpler form: 

W^uv = e^ Vw{W) , Z„, = e^ Vz{W), (38) 

with certain "potentials" IV(W^) and Vz{W). But the simplification does 
not come without a cost: These potentials are explicitly obtained as functions 
of R rather than W. One finds 



Vi 



w 



and 



V. 



K 



-^ -^/^ -AZ{R) 

' ^R{R-K) 
R-K/2 



^/R{R - K) 



,-AZiR) 



(39) 



(40) 



Despite this disadvantage, the form (38) allows an effective treatment of the 



semiclassical dynamics, as will be demonstrated below. 



Notice that the gauge transformation (16) of i? and S carries over to W 
and Z respectively: 



W' = W, Z' = Z-ln^-ln^ 

du dv 



(41) 



The value of Z in a specific w-gauge will be denoted by ^[...], in full analogy 
with our notation for S. 
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3.2 Large-R asymptotic behavior 



From now on we shall consider a macroscopic BH, Mq ^ K, and restrict the 
analysis to spacetime regions where R^ K. We expand Eqs. ( 35p7 ) to first 
order in the small quantity K/R: 



W = R 



-^ — (I 



(42) 



AZ 



K ^(K 
'4i? + ^U 



The inverse function R{W) is given at this order by 



R = W 



^+2W^^^ + K^) 



(43) 



(44) 



The potentials Viy(iy), Vz{W) can now easily be expanded to first order in 
K/W, 



V. 



w 



^-4i^ + ^W 



Vz 



w 



— of— 



(45) 



3.3 The flux-conserving system 

We now proceed with the above large- i? approximation, keeping only terms 
up to first order in K/R (or K/W)-, Thus we analyze the hyperbolic system 



W^uv = e^Vw{W) , Z,, = e^Vz{W) 



with 



where 



Vx 



w 






K 



(46) 
(47) 



Since Vz = dVw/dW, Eqs. ( 46||47 ) constitute a flux- conserving system. 
This concept was first introduced in Ref. |llj and was later described in more 
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detail in Refs. [12] and [TU]. Our system corresponds to F = Vw = — 1 — g/4, 
hence the generating function [15J is 

ho{W) = W + qlnW. (48) 



In conjunction with the system (46||47) we shall also use the transforma- 



tion between (i?, S) and (W, Z) chopped at first order in g, namely 

W = R-2q\nR, (49) 

R = W + 2q\nW, (50) 

^-S-|^S-^, (51) 



We also rewrite here the shell initial conditions ( 30pT ) in the new variables 
W and Z, to first order in q: 

Wo{U) = W{U,Vo) =Ro- 2q\nRo (53) 

and 

zS'\u)^Z[u]{U,Vo)=Vo-^. (54) 

Ingoing Vaidya-like solutions 



Like any fiux-conserving system, the system (46|[47) admits Vaidya-like so- 



lutions — namely, solutions with a single flux [15j. Each Vaidya-like solution 
is endowed with a mass-function — a function of one of the null coordinates 
which encodes the information about the flux. As it turns out, at the lead- 
ing order an evaporating BH may be approximated by an ingoing Vaidya-like 
solution with a linearly-decreasing mass function. [However, a first-order out- 
going component has to be superposed on it in order to precisely match the 
initial conditions, as we discuss below.] Assuming that the BH was created 
by the collapse of a shell of mass Mq, which propagated along the null orbit 
V = Vo, the appropriate mass function takes the form M{v) = Mq — g(f — fo). 
This choice is motivated by the well-known fact [5] that a 2D macroscopic BH 
evaporates at a constant rate q = N/48; and the suitability of (the approxi- 
mate solution derived from) this mass function is verified in Sec. [slbelow. By 
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shifting the origin of v such that Vq = —Mq/q, we obtain the mass function 
in a more compact form: 

M{v) = —qv = rriuiv). (55) 

The restriction to a Vaidya-hke solution leads to a great simplification, 
because the problem now reduces to that of solving an ODE (rather than a 
system of PDEs). Thus, as described in Ref. |10| . W{u, v) is now determined 
from the ODE W^^ = hQ{W) — M{v) (applied along each line u = const), or, 
more explicitly, 

W,^, = {W + q\nW) + qv. (56) 

The other unknown Z{u,v) is then given by 

Z = ln(-W^,„). (57) 

Alternatively [lOj Z may be obtained from the ODE 

Z,. = l + ^. (58) 

As was mentioned above, we set the origin of v such that the collapsing 
shell is placed at vq = —M^/q. Therefore the parameter vq is negative. 
Also, throughout this paper we assume that the BH is macroscopic, namely 
Mq ^ g, and this implies vq ^ —1. 

Weakly-perturbed Vaidya-like solution 

The above mentioned ingoing Vaidya-like solution well approximates many 
aspects of the CGHS spacetime. However, it fails to precisely match the 
initial conditions at the shell. This mismatch is small, oc {q/W), yet it must 
be fixed in order to satisfactorily handle some of the more subtle aspects 
of the solution (most importantly, the Hawking outflux at FNI). Thus, we 
must fix the ingoing Vaidya-like solution by adding to it a small outgoing 
component, seeded by the mismatch at the shell. Nevertheless, owing to the 
small oc {q/W) magnitude of the mismatch, it will be possible to treat this 
outgoing component as a small (linear) perturbation on top of the ingoing 
Vaidya-like solution discussed above (to which we shall refer as the "core 
solution"). 

In the next section we shall proceed with analyzing the ingoing Vaidya- 
like core solution. Then in Sec. |5]we shall construct the perturbing outgoing 
component and thereby complete the construction of the approximate solu- 
tion. 
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4 Constructing the ingoing Vaidya-like core so- 
lution 

4.1 Processing the ODE for W and introducing H 



The term glnVT in the right-hand side of Eq. (56) makes this equation hard 



to analyze. In order to ease the analysis we define the auxiliary variable 
H = W^v + g, namely 



H = W + qhiW -m^ + q. 



(59) 



The inverse function W{H) cannot be expressed in a closed exact form; 
However, restricting the analysis to first order in q we may use the relation 

W = H + m^ - q[\n{H + m^) + I]. (60) 

Differentiating now H using Eq. ( 59 ) and W^^ = H — q we find 



Further substituting Eq. (60) in the right-hand side and omitting all oc q 



terms, we obtain the ODE for H in its more compact form: 



H, 



HI 



q 



H — qv 



(61) 



The initial conditions for H is to be specified at the shell's orbit, the line 
V = Vq (this line serves as a characteristic initial surface for the non-trivial 
piece V > fo of the CGHS spacetime). We denote it Hq{u) = H{u,vo). 
In Sec. |4.4| we calculate Hq and show that it decreases linearly with the 



Kruskal coordinate U. Then in Sec. u^we express Hq(u) in a few other useful 
gauges. Note that the dependence oi H on u only emerges through the initial 
condition Hq{u). Correspondingly we shall often express this parametric 
dependence in the form H{v\ u), and sometimes ignore it altogether and use 
the abbreviated notation H(v), for convenience. 



We are unable to solve the ODE (61) analytically. However, an approxi- 



mate solution is given in Appendix Kl Furthermore, several exact key prop- 



erties of the solutions of Eq. (61) are described below. 
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Note the advantage of the ODE of H over Eq. (56) for W^^, particularly 
at large H. The latter corresponds to future or past null infinity, where, as it 
turns out, H (like W) grows exponentially in v. The ODE then reduces to the 
trivial one Hy = H+q, which is easily solved (leading to the above-mentioned 
exponential growth). In the ODE for W, on the contrary, one faces the In W 
term, which complicates the asymptotic behavior at future/past null infinity. 

4.2 Some exact properties of H{v;u) 



Given the ODE ( |6l[ ), the function H(v;u) is uniquely determined by the 
initial- value function Hq{u) = H{u, Vq) [this function is explicitly constructed 
below; cf. Eqs. (67||68)]. We therefore start our discussion here by mentioning 



two key properties of Hq{u) which are important for the present analysis: 
First, Hq{u) is monotonically decreasing. It is positive at early u but becomes 
negative afterwards. It vanishes at a certain u value which we denote u^°'^ 
(in Kruskal gauge it corresponds to f/ = —qe^'"°). Second, Hq{u) + Mq is 
positive at any u. 



We first note that the ODE (61) admits a trivial solution H(v) = 0. It 
immediately follows that H vanishes along the line u = u^°^ (but nowhere 
else). 

Next, we define (off the line u = u'^°^ where H vanishes) 

l{v;u) = \n\H{v;u)\. (62) 

It satisfies the ODE 

/,. = ! + TFT— = 1 + Tl^ ' (^3) 

H + rriy ±e' — qv 

where the "±" sign reflects the sign of H. 



Both equations (63) and (61) develop a singularity whenever H + rriy 
vanishes, but are regular otherwise. Since Hq + Mq > 0, the quantity H + rrit, 
is always positive at f = f o and its neighborhood. In principle there could be 
two possibilities, which may depend on u [through the initial value Hq{u)]: 
(i) H{v; u) is regular throughout v > Vq, or (ii) H{v; u) becomes singular 
[H + niy vanishes) at a certain finite v = Vsing{u) > fo. n To treat both 
cases in a unified manner we shall say that the solution H{v; u) is regular 



^As will become obvious later, option (i) occurs at u < u^""^ and option (ii) at u > 
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u 



and well-defined throughout the domain Vq < v < Vf{u), where Vf{u) = oo 
in case (i) and Vf{u) = Vsing{u) in case (ii). 

We can now deduce the following exact properties of H{v; u), which hold 
throughout the domain of regularity v < Vf{u): 

(a) As was already mentioned above, H vanishes along the line u = u^°^. 
This holds throughout the range t>o < f < 0. [At {u = u'^°^, f — )■ 0) the ODE 
becomes singular because H + m^ vanishes.] H cannot vanish at any other 
value of M, because otherwise Hq{u) would have to vanish too at that specific 
u 7^ u^""^ (which is not the case). 

(b) Consequently, H{u,v) has the same sign as Hq{u) at any u. As was 
mentioned above, this sign is positive at -u < u^°^ and negative at m > u^°^. 
These two domains are separated by the line u = u^°^ on which u vanishes. 

(c) Since H + m^ is positive at f = fo, it remains positive throughout 
f < ^^ < Vf{u). 



(d) Correspondingly, the right-hand side in Eq. (63) is strictly positive, 
and in fact > 1. Thus, l{v) is monotonically increasing, and the same for 
\H{v)\. 

(e) The quantity Hu satisfies the ODE 

In ii„ = 1 - 



dv ' ' ' {H + m^Y 



[cf. Eq. (87) below]. The right-hand side is regular at any v < Vf{u), and 
positive throughout f < 0. Since i^o.n is negative for all u, f\ it immediately 
follows from this ODE that H^ is everywhere negative. Furthermore, at least 
at f < 0, |-ff,u| is an increasing function of v; therefore, at fixed u (and f < 0) 
Hu is bounded above by the parameter Hqu < 0. 

(f) The same obviously applies to {H + m^)^u (in particular, this quantity 
is everywhere negative). It then follows that if a singularity H + m^ = 
occurs at some u = Ui (with finite Vsing > Vq), then such a singularity must 
occur at any u > Ui, and Vsmgiu) must be a non- increasing function of u 
(throughout u > Ui). Furthermore, in the range where Vsing{u) < [see (g)], 
Vsingiu) must be a strictly-decreasing function of u. 

(g) As was mentioned above, at m > u^°^ Hq{u) is negative, and the 
same for H{v;u). [Furthermore, since |-ff(f)| is monotonically increasing, at 
each line u = const in this domain H is bounded above by the parameter 

^Hq decreases linearly with U; and we only consider here u-gauges satisfying du/dU > 
0. 
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Hq(u) < 0.] It then follows that H + m^ (which is positive at v = Vq) must 
vanish before rriy vanishes. Thus, all lines u > u^""^ run into an if + m^ = 
singularity at a certain v = Vsing{u) < 0. Property (f) then implies that 
Vsing is a strictly-decreasing function of u, indicating that this is a spacelike 
singularity. 

(h) In the domain u < u^°^, H is strictly positive, and no singularity 
may form at t; < (where rriy is positive too). However, at f > m^ is 
negative, and one might be concerned about the possibility of vanishing H + 
rriy, which would lead to a singularity. A closer examination reveals that such 
a singularity does not occur in this domain. This may be deduced from each 
of the following arguments (though a complete mathematical proof is still 
lacking): (i) It is easy to show that at least a locally-monotonic singularity 
of this type (namely, a singularity of vanishing H+m^ at some finite v = Vging, 
such that H + m^ is monotonic in v throughout some neighborhood oi v = 
Vsing) IS not possiblc in the domain u < u^"^: H+rUy starts positive at f = Vq, 
and in order to vanish it must decrease on approaching v = Vging- However, 
since H > (and increasing), when H + ruy approaches zero from above Eq. 



(61) yields H^y — )■ +00, and therefore H + rriy must increase, so it cannot 
vanish. It is harder to mathematically exclude the possibility of an oscillatory 
approach to an ii^ + m^ = singularity, but such an oscillatory behavior 
seems very unlikely, (ii) Consider the limiting function H{v;u'^^), defined 
to be the limit u — )■ u'^'^ of H{v;u). At f < this function vanishes, just 
like H{v] u^°^), by continuity. At v > [where H{v]u^°''') is not defined; see 



(a)], this function becomes a non-trivial solution of the ODE (61 ). Numerical 
examination shows that H{v]u^'^) continuously increases from zero at f = 
to infinity at f — )■ 00, with H + my > at any f > 0. From (f) it now follows 
that H + rriy > throughout u < u^°^. (iii) Direct numerical simulations of 
H{v;u) at various u < u^°^ values further confirm this conclusion. 

Summarizing the above discussion on the exact properties of H{u, v), and 
briefiy re-stating it in more physical/geometrical terms: All lines u < u^°^ 
make it to FNI, whereas all lines u > u^°^ run into a spacelike singularity |j 



^The exact semiclassical CGHS spacetime admits a very similar global structure: A 
BH, with a spacelike singularity inside it. We point out. however, that the local properties 
of the spacelike singularity in our approximate solution are different from those of the 
precise CGHS spacetime. In particular, in our approximate solution (when literally applied 
to the H + m.„ = singularity) R diverges logarithmically to —00, whereas in the CGHS 
solution R = K 'At the singularity |7j. This remarkable difference is not a surprise, because 
the parameter q/R is no longer small when we get close to the singularity, hence our 
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&tV = Vsingiu) < 0. 

Based on these causal properties of H{u, v) and its singularity, we shall 
refer to the ranges u < u^°^ and u > u^°''' as the BH exterior and inte- 
rior, respectively. Note that the interior is confined to the range Vq < v < 
[u) < (and u > u^"^), whereas the exterior extends in the entire domain 



-'sing 



■^0 < "^ < oo (for u < u 



hor\ 



4.3 Expressing the ingoing solution in terms of H 



The function H{u, v) may be obtained by solving the ODE (61 ) numerically, 

or by analytic approximate solutions like the one given in Appendix [Xj Once 

H is known, W is given by Eq. (60), and Z in turn by Eq. (57). The former 

equation yields 

Q 



W,u = (1 
hence, up to first order in q, 

Z = H-H,u 



H + TUy' 



)H, 



q 



(64) 



H + niy 
One can easily verify the consistency of this expression with the w-gauge 



transformation of Z, given in Eq. (41) (note that H is unchanged in such a 
transformation) . 



4.4 Calculating Hq{u) 



To calculate Hq{u) we evaluate Eq. (64) at v = vq and equate it to the 



desired initial condition Zq{u). It is convenient to carry out this calculation 
in the Kruskal gauge. Using Eqs. ( 54p0 ) we obtain the following equation 
for Ho: 

^"^-^°'^) - luUio = "° - M^^- ('') 

A simple solution immediately suggests itself: Hq = — f/e^"; but we need here 
the general solution for this ODE, which spans a one-parameter family. For 
our purpose it will be sufficient, however, to derive an approximate general 
solution (up to order q), which is an easy task. Such a one-parameter family 
of approximate solutions is 



HoiU) = -Ue'"'+pq, 



(66) 



approximate solution becomes invalid there. 
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where p is a yet-arbitrary constant. This constant may be fixed by solving the 



ODE (61 ) for H{v; U) in the PNI asymptotic hmit [with the above expression 



for Hq{U) as initial condition], constructing W, Z from H and then i?, S", and 
comparing them to the desired initial data at PNI. In Appendix [B] we carry 
out this analysis and find that the appropriate value is p = — 1 , namely 



Ho{U) 



-Ue 



VQ 



Q- 



(67) 



We may also use Eq. (|30|) to re-write Hq in a form which is explicitly gauge- 

(68) 



invariant (for arbitrary u coordinate): 

Ho{u) = Ro{u) -Mo-q. 



5 The weak perturbing outflux 

The above ingoing Vaidya-like solution was constructed (through an appro- 
priate choice of Hq{u)) such that it properly matches the initial function 
Zo[u) at the shell. However, the other initial function Wo{u) does not ex- 
actly coincide with W of the above constructed ingoing solution. To fix this 
mismatch we shall add a weak, outgoing component as a perturbation on top 
of the ingoing core solution. 

To this end we write the overall W, Z functions as 



W{u,v) = W'''{u,v) + 6W{u,v), (69) 

Z{u,v) = Z'''{u,v) + 6Z{u,v), (70) 

where W^"', Z*" denote the Vaidya-like ingoing core solution constructed in 



the previous section [namely Eqs. (60) and (64)], and 5W,6Z denote the 
additional perturbing component, q] 

We first calculate the mismatch in the initial condition for W. Equations 



(|60j) and (|67|) yield for VT™ at the shell 

W^'^iU) = Mo - Ue""" - g[ln(Mo - Ue"") + 2] 



(71) 



(we have omitted the q in the log argument, being a higher-order term; and 
we do this occasionally in the equations below). This is to be compared to 



the actual initial data for W, Eq. (53), which, together with Eq. (30), reads 
Wo{U) = Mo - Ue''° - 2q ln(Mo - t/e'^o). (72) 
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Note that in a coordinate transformation Z™ transforms like Z, hence SZ is invariant. 
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The difference is thus 

6Wo{U) = -g[ln(Mo - f/e^°) - 2]. (73) 

Note that no mismatch is present in the shell data for Z, because we have 
chosen Hq{U) in the first place so as to properly match Zq{U); therefore, 

SZoiU) = 0. (74) 

Next we substitute Eqs. (69|[70) in the full flux-conserving system (46), to 



obtain field equations for 6W, 6Z. Since we are only interested in the solution 



up to first order in g, and the mismatch initial data (73) are already 0(g), we 
may treat 6W, SZ as small perturbations, satisfying the linearized equations 

6W,uv = e^ [Vw,wSW + VwSZ] , (75) 

6Z,,v = e^ [Vz,wSW + Vz6Z] . (76) 

Furthermore, we only need to consider here the coefficients (Vw, Vw,w, Vz, Vz,w) 
at zero order in q — namely, Vw = — 1 and Vw,w = Vz = Vz,w = 0, yielding 
the trivial system 

SW^uv = -eHZ , bZ^^, = 0. (77) 

The initial conditions are Eqs. (73174) at the shell, and no contribution from 



PNI. P^For 5Z we immediately obtain 

6Z{u,v)=0. (78) 

In turn 6W satisfies the trivial equation 6W^uv = 0, yielding 
6W{U,v) = 6Wo{U) = -g[ln(Mo - f/e^") - 2]. 
Thus, the overall solution is 

W{U,v) = W'''{U,v) + 6Wo{U) = i^ + m„-g[ln(/7 + mJ + ln(Mo-f/e^")-l] 

(79) 
and 

Ziu, v) = Z'^'iu, v) = H-H^) - jr^- (80) 

SI + 777,,, 



Finally we re-write Eq. (79) in a form which includes no specific reference 



to the Kruskal coordinate U, by replacing Mq — Ue^° (in the log argument) 



with Ho + Mo, using Eq. (67): 



W{u, v) = H + m,- q[\n{H + m,) + \n{Ho + Mq) - 1]. (81) 



^^ Since the perturbation we add here is an outgoing component, it is assumed to be 
seeded at the shell only (any ingoing component would be absorbed in the ingoing core 
solution in the first place). 
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6 The final approximate solution 



Transforming the above results ( [SIpOJ ) from W, Z to the original variables 
R, S, using Eqs. (50152), we obtain our approximate solution in its final form: 

(82) 



R{U, v) = H + m, + q[\n{H + m,) - \n{Ho + Mq) + 1] 



and 



S{u,v) = \n{-H^^), (83) 

where q = K/A. The function H{u^ f), recall, is determined by the ODEP] 

^ ' (84) 



(85) 



\ H — qv 
with initial conditions Hq{u) = H{u,Vo) given by 

Ho{u) = Ro{u) -Mo-q 
or, more explicitly (in the Kruskal gauge) 



HoiU) 



-Ue 



Vo 



(86) 



Approximate analytic expressions for H{u,v) are given in Appendix^ cf. 
Eqs. p7p08| ). 

Note that this approximate solution precisely matches the required initial 
conditions at the shell, namely (using Kruskal [/-gauge) R = Mq — Ue'"° and 
S[u] = \n{-H^u) = Vo- 



The expression (83) for 5* requires Hu- The approximate analytic expres- 
sions ( 107|108 ) can be directly differentiated to yield H u- However, when 
H{v; u) is obtained by numerically solving the ODE (84), a direct numerical 
^-differentiation may be inconvenient. In this case it is easier to obtain H^ 
by numerically integrating the ODE it satisfies: 



dv 



{H,u) 



1- 



2 

q V 



{H-qvy 



H, 



(87) 



(this can be done simultaneously with the numerical integration of the ODE 
of H itself). The initial condition at the shell is obviously 



HJu,vo) 



d 
du 



Ho{u). 



(88) 



^^Note that the parameter q may easily be factored out of this ODE: Defining H = H/q, 
we obtain the ODE in its universal form H,v = H[l + 1/{H — v)]. 
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6.1 Alternative approximate expression for S 

As was already mentioned in Sec. [3} the Z function in the ingoing Vaidya-hke 
solution can be obtained by either of the equations (57) or (58). The above 



analysis was based on the former equation, and it led to the expression (83) 
If one uses Eq. (58) instead, one can derive an alternative expression for S: 



galt^ 



U,V 



Vq + ln 



H 



Hoiu] 



+ q 



1 



1 



H + m^ Ho{u) + Mo_ 



In 



du 
dU' 



(89) 



This expression for S looks more complicated, but it has the advantage that 
it does not require H^- 

It should be emphasized that the two expressions ( 83p9 ) for S are not 
exactly identical. Yet the difference appears to be compatible with the an- 
ticipated error characterizing the entire approximation scheme used here. At 
the same time we also point out that so far the error in the approxima- 



tion (89) for 5* has not been explored as thoroughly as that in the original 
approximation (83) (cf. Sec. [s]). 



7 Some useful gauges 



Our approximate solution was presented in Eqs. (82 85) in a fully (M-)gauge 



covariant form. The only reference to a specific gauge was made in Eq. (86), 
which explicitly gave Hq in terms of the Kruskal [/-coordinate. 

In this section we shall introduce a few additional useful w-gauges: The 
shifted-Kruskal gauge, and the (external as well as internal) semiclassical 
Eddington gauge. These new gauges slightly simplify the functional form of 
Hq{u). More importantly, they are better adopted to the global structure of 
the evaporating-BH spacetime (e.g. the location of the horizon). 

Note that in all these gauges, we use for the f-gauge the same Eddington 
coordinate v (originally denoted Vg), as we do throughout this paper. 



7.1 Shifted-Kruskal gauge (U) 

We define the shifted Kruskal coordinate 

U = U + qe-''\ 
Since this is a constant shift, S is unchanged: S'rjji = S[u]. 
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(90) 



The expression for Hq slightly simplifies in this gauge: 

Hq{U) = -Ue"\ (91) 

Note that U^°^ = 0. 

We point out that the solution in {U, v) coordinates [just like in {U, v) 
coordinates] covers the entire BH spacetime. As may be obvious from the 



discussion in subsection 4^, the BH exterior and interior correspond toll < 
and U > 0, respectively. 

7.2 Semiclassical Eddington gauge (u) 

In the range U < (the BH exterior) we define the Semiclassical Eddington 
coordinate u by 

u = -ln{-U) {U<0). (92) 

S is modified by this transformation according to S^u] = S^u] — u. The initial 
function for H now reads 

Hq{u) = e"°-^. (93) 



In the Semiclassical Eddington gauge the alternative expression (p9j) for 
S reduces to the simpler form 



Sl^{u,v)=\nH + q 



1 1 



(94) 



H + m^ Hq + Mq 
Note that the coordinate u only covers the BH exterior. 

7.3 Internal semiclassical Eddington gauge {u) 

In the range f/ > (the BH interior) we define the internal semiclassical 
Eddington coordinate u by 

u = \n{U) {U > 0). (95) 

Now S is modified according to S^u] = S[u] + u. The initial function for H 
takes the form Hq{u) = e^°~^^. 
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8 Verification of the approximate solution 



Our approximate solution (82 85) was constructed here through a rather 
indirect process, which involved the transformation to new field variables 
W, Z, the large- i? approximation, the formalism of flux-conserving systems, 
and their Vaidya-like solutions. It is therefore important to directly examine 
the validity of the resultant expressions. 

Naturally this examination involves two independent parts: (i) checking 
compliance with the field equations; and (ii) checking compatibility with 
initial conditions, both at the shell (f = Vq) and at PNI. These two parts will 
be carried out in the next two subsections. 

8.1 Compliance with the field equations (error estimate) 

To define the local error in the evolution equations we substitute the ap- 
proximate expressions (82p3) in the field equations (T3p4) and evaluate the 



error — namely, the deviation of R^uv and S^uv from their respective values 
(specified at the right-hand side of these two equations). The error defined 
in this way is obviously gauge-dependent, and we find it convenient to employ 
the semiclassical Eddington coordinate u (along with Eddington v) for this 
task. Using the MATHEMATICA software we find that the local errors in the 
two equations indeed scale as g^, as anticipated. More specifically, the errors 
scale as R{q/ RY for R^^y and as [q/ RY for S^uv — both multiplied by certain 
functions of rriy/R. This local error estimate applies to typical off- horizon 
strong- field regions, namely, regions for which m^/R is of order unity but not 
too close to 1. The error decays exponentially in 1^; — «! both at weak-field 
regions {m^/ R — )■ 0, corresponding to f — -S — t- oo) and near-horizon regions 
{rriy/R — )■ 1, corresponding to f — u — )■ — cxd). 

The global, accumulated, long-term error is harder to analyze. It may 
be evaluated by comparing the approximate solution (82p3) to numerical 



simulations, but this is beyond the scope of the present paper. In the next 
section, however, we shall evaluate the accumulated error in both R and 5* 
in the neighborhood of the horizon (m — )■ oo). 

The error in the constraint equations (9p0) is also found to be propor- 



tional to g^, as may be expected (based on the mutual consistency of the 
evolution and constraint equations). However, the functional dependence of 
the pre-factor on m^/R is more subtle and will not be addressed here. 
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8.2 Compatibility with the initial conditions 
Initial data at the shell 



At f = f the expressions (82) and (83) reduce to R = Hq + Mq + q and 



S = \n{—HQu) respectively. By virtue of Eq. (86), one obtains (using the 



Kruskal U coordinate) R{U,vo) = Mq — Ue"° and S[u]{U^vq) = vq, which 
exactly match the desired initial data (30pl). 
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Initial data at past null infinity 



In Sec. 10 we analyze the asymptotic behavior of our approximate R and 5* 



at PNI, and verify that they satisfy the required asymptotic behavior (32). 



We also show that the influx T„^ = RyS y — R^y vanishes at PNI, as it should. 



9 Horizon 



In Sec. H^ we analyzed the behavior of H{u, v) along u = const lines, and 
found that spacetime is divided into two domains by a certain outgoing null 
ray which (for a general m- gauge) we denoted u = u^°^: All lines u < u^°^ 
run to FNl in a regular manner (with steadily growing H), whereas all lines 
u > u^°^ crush into a spacelike singularity. The spacetime thus contains a 
BH, and the domains u < u^°^' and u > u^""^ correspond to the BH exterior 
and interior, respectively. We shall therefore regard the critical null ray 
u = M^°'' as the event horizon (or sometimes just horizon) of the BH . ^^ 

The horizon is characterized by the vanishing of H, hence its location 
u = u^°'^ is determined by requiring Hq{u) = 0. Referring to some specific 
gauges, the horizon's location is [/ = in the shifted Kruskal gauge, and 
U = —qe~'"° in the original Kruskal gauge. In the semiclassical Eddington 
gauge the horizon is located at the asymptotic boundary m — )■ oo. 

Recall that in the classical solution (with the same initial data) the hori- 
zon is located at U = 0. The inclusion of semiclassical effects thus shifts the 
horizon in U (or U) by an amount qe~'"° . 



^•^Obviously, this also implies exact matching of R and S to the Minkowski-Hke solution 
at w < uq in any other gauge. 

^''This involves some abuse of standard terminology, because the line u = u^""^ actually 
contains a naked singularity at w = (the point denoted "P" in Fig. 1), where H + m„ 
vanishes. Note that it is only the section v < of this line which separates the BH interior 
and exterior. The portion v > Q oi u = u is in fact a Cauchy horizon (see Fig. 1). 
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9.1 Behavior of R and S at the horizon 

The behavior of R along the horizon is obtained by setting Hq = H = in 



Eq. (82), yielding rriv + g[ln(mi,/Mo) + 1]. However, the term q\n{mv/Mo) 
in this expression cannot be trusted, as may be deduced from simple error 
estimate. To this end we evaluate the error in R.^ as a function of v along 
the horizon, using semiclassical Eddington coordinates for simplicity. We 
do this by integrating the local error in Ru^ along a line v = const < 0, 
from u = — oo (PNI) to tt = oo (horizon). From the discussion in Sec. M 
it follows that along such a v = const line the local error gets a maximal 
value of order q'^/Rr^ q'^/ni^ at intermediate u ~ —v values, and it decays 
exponentially in u in both directions. The effective integration interval (in 
u) is of order unity, hence the integrated error in R^^ is also ~ q'^/niy = q/\v\. 
Consequently, the integrated error in R{v) along the horizon is of order ~ 
gln(f/fo) = gln(m^/Mo). rj We therefore re-write the above result for R{v) 
asE3 



R^'^'iv) =m„ + g + Ogln— ^ = -qv + q + O [qhi— ] . (96) 

Next we analyze the behavior of 5* in the horizon's neighborhood, in 



the semiclassical Eddington gauge, using Eq. (83) which now reads S\.i^ 



ln(— if ii). To this end we divide Eq. (87) by H u, substitute if = in the 



right-hand side, and re-write this equation as 

(97) 



d , 




1 


, In Eu)- 


1 - 




dv ^ ' ' 




V 



The initial condition at f = f o is obtained from Eq. (93) which yields 



ln(— iio.u) = Vq — u. Integrating Eq. (97) we find 

b\u] = f — M -|- in — . 

V 

The error estimate for S parallels the one carried out above for R. It yields 
an accumulated error ~ {q/rriy)^ = v""^ in the value of S^^ at the horizon. 



^^It should be emphasized that despite this 0{q) integrated error, it is crucial to keep 
the 0{q) term in the approximate expression ( [82| for R. Without this term, the local 
error in R^^ will grow from 0{q^) to 0{q). 

^^In fact one can use the constraint equation for r^„ to obtain the correct coefficient of 



the log term in Eq. (96 I, but this is beyond the scope of the present paper. 
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and hence an integrated error of order ~ ^/\v\ = q/m^ in S itself. W\ We 
therefore re-write our resuh as 



-hor _ „, „-, , . ^'0 , ^,, ,. _ „, ~ , , Mo 



Sf^ = v-u + \n— + 0{l/v) = v~u + \n^^ + 0{q/m^). (98) 

' ' V ruv 



Another way to obtain this resuh is by integrating Eq. (63) for /(f). In 
the horizon's neighborhood this equation reduces to l^^ = 1 — 1/v, which 
is easily integrated (with the appropriate initial conditions) to yield / = 

V — u + ln(fo/f), or 

ffhor ^ ^^v-u_ /99^ 

V 

Differentiating this expression with respect to u and substituting in Eq. (ISSJ), 



one recovers Eq. (98). 
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9.2 The apparent horizon 

Motivated by the terminology used for conventional 4D spherically-symmetric 
BHs, we define the apparent horizon to be the locus of the points where 



R^v = (or equivalently, 0^^ = 0). By virtue of Eq. (82) this implies 



(H + m^)^v = 0, or Hy = q. Utilizing Eq. (84), and restricting the analysis 
to first order in q, we find that at the apparent horizon H = q is satisfied. 

We shall consider here the properties of the apparent horizon during the 
macroscopic phase m„ 3> q, namely v <^ —1. Throughout this phase H^°^ ^ 
R = rriv and hence the horizon approximation H = H^"^ applies. Setting 



if = g in Eq. (99) we find the apparent-horizon's location 

u^v + ln{vo/v)-\nq. (100) 

Thus, du/dv = 1 — 1/v ~ 1 along the apparent horizon. We find that the 
apparent horizon is a timelike line, which is approximately vertical in the 
(u, f ) coordinates. In Fig. 1 the apparent horizon is denoted "AH". 

Consider now the behavior of H and R along an outgoing null geodesic 
located outside the BH though fairly close to the event horizon — namely, 
sufficiently-large fixed u. To be more specific, let us assume that fo < m < 0, 
such that -u ^ — 1 but u — vq^ 1. (As a typical example one may take u ~ 



^^To be more precise, the integrated error is ^ {'^/\v\ — l/|wo|) — [ql'm-^ — q/Mo). 
^^The alternative approximation S""'* yields the same result. To see this one substitutes 



Eq. (99 1 in (94 1 (and omit the unsecured 0{q) term). 

31 



fo/2; Recall that throughout this paper we assume Vq ^ —1, corresponding 
to a macroscopic BH.) Then Hq = e^°~" is exponentially small and may be 



neglected in Eq. (82). Initially, near v = Vq, H is also exponentially small 
and therefore (neglecting terms of order ~ q compared to m^,), R ~ m„. In 
this range R shrinks linearly in v, like m„ itself (and like R'^"^). During this 
stage H = H^°^ grows (approximately) exponentially in v, but initially this 
does not have much effect on R{v) because H is still too small. However 
at some point the exponentially-growing H starts to slow the decrease rate 
of R{v). When H approaches q this exponential growth just balances the 
linear decrease of rriy. This is the point of intersection with the apparent 
horizon. Note that at this point R is still ^ m^, the difference being 0{q). 
Soon afterwards the exponentially-growing H overtakes m^,. 

On the other hand, for earlier outgoing geodesies with sufficiently large 
Hq, the exponential growth of H will dominate over rra^, everywhere, and R 
will grow monotonically all the way from vq to oo. Since H,^ > outside the 
BH, and the apparent horizon satisfies H = q, this monotonic growth of R 
will occur at the outgoing null geodesies for which Hq > q. 

The evaporating-BH spacetime may thus be divided into three domains 
in u. For concreteness let us use here the coordinate U to characterize these 
domains: In the early domain, U < Uah, R grows monotonically with v 



throughout Vq < v < oo. From Eq. (91) we find (equating Hq to q as 
explained above) 

Uah = -qe-""'. (101) 

In the second domain, Uah < U < 0, R first decreases along an outgoing 
null ray (taking values fairly close to my) until it intersects the apparent 
horizon at certain v < 0, and then R starts to increase with v. In the third 
domain U > (the BH interior), H is everywhere decreasing (and the same 
for rriy), therefore R decreases monotonically until the outgoing null ray hits 
the spacelike singularity. 

From the discussion above it follows that at a given f < the event 
and apparent horizons have roughly the same R value (the difference being 
approximately q). In other words, the apparent and event horizons shrink 
(in R) in the same standard rate, dR/dv = —q. 

It may be interesting to compare three points along the worldline of the 
collapsing shell: (1) The intersection oi v = Vo with the apparent horizon 
(f/ = If ah)] (2) its intersection with the event horizon {If = 0); and (3) 
its intersection with the (would-be) event horizon of the classical CGHS 
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BH (namely U = 0, or U = ge~^°). One finds that these three points are 
equally-separated in U (or in U), the separation being qe~^°. We find it 
more illuminating, however, to express these three points by their respective 
R = Ro{u) values. Noting that 

one finds that Rq is Mq + 2q at point 1, Mq + g at point 2, and Mq in point 
3, so these points are equally-separated in R too. Notice, however, that 
for a macroscopic BH (Mq ^ q, which we assume throughout) the relative 
separation g/-Ri,2,3 ~ q/Mq is ^ 1. 

9.3 Influx at the horizon 

We proceed now to calculate T^„ at the horizon, using the constraint equation 



(10) which now reads 

rpnor Tjtior qtior r), 



hor T^hor nhor r)hor 

vv \v ,v \vv ' 



We shall restrict here the calculation to the leading order, namely first order 
in q/m^. From Eqs. (|96|98|) we find S'^^"'' = 1 + 0{q/m^), R^^°^ = -q + 



0{q'^/my), and R'^°^ = U[q /ml) will not contribute. We obtain at the 
leading order 

r,r = -9 = -f = -^. (102) 

Thus, in the macroscopic limit the infiux into the horizon is constant and 
independent of the mass — a well-known result for a two-dimensional BH |2j . 



In principle it is possible to use Eq. (12) for T^j, to obtain the first-order 



correction to the fixed infiux (102) (anc. 



thereby to fix the O[gln(m„/Mo)] 



term in R^°^), but this is beyond the scope of the present paper. 



10 Null infinity 

10.1 Past null inflnity 

The limit U — t- — oo (also U,u ^ — oo) and finite v corresponds to PNI. 
In this asymptotic boundary m^ is finite but H diverges (this immediately 
follows from the divergence oi Hq oc —U, combined with the monotonic 
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growth of |ilf| with v) 

H + q and its solution, corresponding to the initial conditions (|86|) 



At this limit the ODE ^^ for H reduces to H^^ 

is 



HPm ^ _f/g^ 



Q- 



The expression (83) for 5* then yields 



S 



pm r^ 



[U] 



(For the other gauges we find S^^^ = v and 5*?^^ 
yields R 



pni 



U. 



=m„ 



Ue^ + q\n{H/Ho), and setting \n{H / Hq) 



(103) 



(104) 

For R Eq. (|82|) 
f — I'd we obtain 

(105) 



These expressions for i?^™ and S*^"* properly match the desired initial con- 



ditions (32) at PNI. 



By assumption the influx T^^ should vanish at PNI. We would like to 
verify this by applying the constraint equation T„^ = R^^S^v ~ R,vv to our 
approximate solution and taking the limit U — t- — oo. Doing so we observe 
that Eqs. ( |T04p05| ) yield R^f = Rp^^ = -Ue" and S^^' = 1, which yields 
the desired result 

TP™ = 0. (106) 

Note the following subtlety, however: Because i?^"* oc U diverges at PNI, 
the term R^^^S^^^ might receive a nonvanishing contribution from 0(1/ U) 
corrections to S, if such existed. To address this issue we must carry the 
calculation of Sf^! to order 1/U. This in turn requires a more detailed exam- 



ination of HP''\ At PNI the ODE 1^ takes the form H^^ = H + q + 0{l/H), 
and the last term introduces 0{1/U) corrections to H, namely ifP"* = 
—Ue"" — q + 0{1/U). How ever, this only leads to 0(l/?7^) corrections in 
Sf^T, which leave Eq. (106) intact. 



10.2 Future null infinity 

The limit f — )■ oo (for U < 0) corresponds to FNI. The analysis of this asymp- 
totic region is more complicated than that of PNI, because this time we have 



to integrate the ODE (84) through strong- field regions (where H is compa- 
rable to rriv). The thorough investigation of this asymptotic region is beyond 
the scope of the present paper, and we hope to address it elsewhere. Here we 
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shall merely mention two key properties: (1) Spacetime is asymptotically- 
flat; In other words, for appropriate choice of u-coordinate, which we denote 
u (and which does not exactly coincide with u), p[^] vanishes at f — )■ cxd. 
(2) At the leading order one obtains a constant, mass-independent. Hawking 
outflux T^^ = q = N/4:8 (a well-known result [2j). 

11 Summary 

Our approximate solution for the semiclassical variables R and S is described 
in Eqs. (82p3). The determination of the original CGHS variables (p, p from 



R and 5* is straightforward. The function H involved in this solution is 



determined by the ODE (84), along with the initial conditions (85) or (86). 



Approximate expressions for H are given in Eqs. (107[108). 



In the CGHS formalism all semiclassical effects originate from a term in 
the action which is proportional to K = N/12. Our approximation scheme 
is restricted to macroscopic BHs, namely, those with original mass Mq 3> 
q = K/4. Furthermore, it only applies to spacetime regions where R ^ q. 
It nevertheless holds both outside and inside the BH, though obviously not 
too close to the singularity (where the condition R^ q is violated). 

Throughout this paper we use double-null coordinates {u,v). Our con- 
struction explicitly preserves the gauge-freedom in u, though not in v. Our 
coordinate v coincides with the affine parameter along PNI. (We have set 
the origin of v such that the collapsing shell is placed at v = vq = —Mo/q.) 
For the w-gauge we flnd two particularly useful choices: the shifted-Kruskal 
coordinate If, and the semiclassical Eddington coordinate u (both deflned in 
Sec. Q. 

At the horizon's neighborhood we flnd that R shrinks linearly with v, 
as may be anticipated: R = —qv + q. The other variable 5* behaves as 
S = V + ln(fo/f ) + c, where c is w-independent though it depends upon the 
M-gauge being used. (For example, c = in the shifted-Kruskal gauge and 
c = —u in the semiclassical Eddington gauge.) The term v merely reflects the 
classical behavior of S at the horizon, whereas the log term is of semiclassical 
origin. 

Our approximate solution is "locally flrst-order accurate"; Namely, the 
deviation of R^^^ and S^u^ from their respective values, as specifled in the 
evolution equations (6|7) or (13p4), is proportional to g^. [More speciflcally. 



the deviations in both R^uv/R and S^ny are ~ {<i/RY at typical (off-horizon) 
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strong-field regions, and smaller elsewhere.] On the other hand, the accu- 
mulated error in S, and accumulated relative error in R, are of first order 
in q/R. This change in the power of q is because the effective accumulation 
intervals (in v and/or u) are typically of order of the evaporation time, which 
is Mo/q. 

It will be interesting to check this approximate solution against numer- 
ical simulations of the CGHS field equations. Such numerical simulations 
are currently being conducted by several groups O [6]. A preliminary com- 
parison with the numerical results [5j shows nice agreement, though a more 
comprehensive check still needs to be carried out. 

ACKNOWLEDGEMENTS 

I would like to thank Liora Dori, Fethi Ramazanoglu and Frans Pretorius for 
helpful discussions, and for sharing their numerical results. I am especially 
grateful to Abhay Ashtekar for his warm hospitality and for many interesting 
and fruitful conversations. 



A Approximate solution for H{v]u 



A useful approximate solution of the ODE (84) for H is 



Hiv;U)^-Ue^'''^'-''-^/'\ (107) 

^ ^ ln(e- - U/q) ^ ' 

where U = U + qe~'"° is the shifted Kruskal coordinate defined in Sec. [tI This 
expression holds both inside and outside the BH (corresponding to t/ > 
and f/ < respectively). 

Note that the initial condition Hq{U) = — f/e^" at f = fo is precisely 
satisfied. Also recall that the dependence oi H on U only emerges through 
this initial condition at the shell. It is thus possible to substitute U = 



—Hoe~'"° in Eq. (107), to obtain an expression for H in terms of Hq{u) but 
with no explicit reference to any specific w-gauge. 

We also translate the above expression for H to the semiclassical Edding- 
ton coordinate -u, valid in the external world ([/ < 0): 

^ ' ln(e-^ + e-"/g) ^ ' 
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Preliminary error estimate suggests the following behavior of the relative 
error in the above expression for H\ We assume fo <^ — 1 throughout. For 
M ~ fo or fo < u ^ — 1 the error scales as 1/m, though with some logarithmic 
corrections. This applies to both (i) the relative error in H itself, and (ii) the 



relative local error in if „ [compared to its respective value in the ODE (84) 
For u < Vq the error further scales as e"~^° and quickly becomes negligible. 



B Calculating the constant p 



In this Appendix we determine the constant p in Eq. (66) by analyzing the 



asymptotic behavior of R at PNI and comparing it to the desired initial 
conditions. 

PNI is characterized by finite v but [/ — )■ — cxd, yielding finite niy but 
H ^- oo (the latter follows from the divergence of Hq combined with the 



monotonic growth of \H\ with v). Therefore, in the ODE (84) we may replace 
the term H — qv in the denominator by H . We are left with the ODE 
H y = H + q, whose general solution is 



H{u,v) = -q+[Ho{u) + q]e^- 



■VQ 



Using the general expression (66) for Hq{U) we get 



H{U, v) = l-Ue^'' + {p + l)g]e''-^« - q. 



(109) 



We now substitute this in Eq. (60 ) for W (or more precisely, W^""; see below). 
In doing so, we may approximate qln{H + m^) at PNI by gin if, which by 



virtueofEq. (109) may be approximated by g[f+ln(—f/)] = — m^+gln(— [/), 
omitting 0{q^) contributions. (The same will apply to terms like glniy and 
glni?o below.) We obtain 



W^'^iU^v) = [-f/e^o + {p + l)q]e''-'"' + 2m„ - 2g - q\n{-U). 



;iio) 



We used here the symbol H^*" to make it clear that this is the W^-function 
associated with the ingoing Vaidya-like core solution (as opposed to the fuU- 
W function, considered below; See discussion in Sec. [s]). At f = f o this 
yields 

W'''{U,Vo) = -[/e^o + 2Mo - gln(-f/) + {p - l)q. 
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This should be compared to the full-VT initial data at f = fo, given by Eqs. 
( |53p0| : 



Wo{U) = Mo - Ue"" - 2q[vQ + ln(-f/)] = 3Mo - f/e^" - 2gln(-f/). 



[Here, again, in processing the term 2qlnRo = 2gln(Mo — Ue^°) in Eq. (53) 
we have neglected Mq compared to Ue""".] The difference is 

6W{U) = Wo{U) ~ W'^'iU^vo) = Mo - q\n{-U) - {p ~ l)q. (Ill) 

The full function W{U, v) is obtained by simply adding the outgoing 
perturbation 6W{U) to the ingoing solution W^^{U, v) (this is shown in detail 
in Sec. [Sl based on linear perturbation analysis). Summing Eqs. ( 110|111 ) 
we obtain for the full-W^ function near PNI 

W = [-Ue"° + {p+ l)g]e^-^" + Mq + 2m^ - 2g ln(-t/) - (p + l)g. 



Finally we transform from W to i?, using Eq. (50 ). The term 2g In W therein 
[in which we approximate In VT = v +\n{—U) as explained above] just cancels 
the terms 2mjj — 2q\n{—U) in the last expression for W, and we obtain 

R = Mo-Ue'' + {p + l)g(e''-^'' - 1). 



This expression should agree with the presumed PNI initial conditions (32), 
namely R = Mq — Ue^, which dictates p = —1. 
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Figure 1: Penrose diagram describing the formation and evaporation of a 
two-dimensional CGHS black hole. The black hole forms by the collapse 
of a thin shell of macroscopic mass Mq. The collapsing shell is located at 
V = Vq = —Mq/q. The lines denoted "EH", "AH", and "CH" respectively 
represent the event horizon, apparent horizon, and Cauchy horizon. 
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